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Abstract 

The dynamics of a one-dimensional crystalline interface model with long-range 
interactions is investigated. In the absence of randomness, the linear response 
mobility decreases to zero when the temperature approaches the roughening 
transition from above, in contrast to a finite jump at the critical point in the 
Kosterlitz-Thouless (KT) transition. In the presence of substrate disorder, there 
exists a phase transition into a low-temperature pinning phase with a contin- 
uously varying dynamic exponent z > 1. The expressions for the non-linear 
response mobility of a crystalline interface in both cases are also derived. 
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Various deposition growth models of interfaces have been investigated intensively in recent 
years |l],|2|. The main interest of these studies stems from both practical applications and 
theoretical inquires. The profile of a surface develops roughness gradually as the stochastically 
deposited particles accumulate. The primary features of interface growth are their non-trivial 
temporal and spatial scaling behavior. The scaling properties are manifested in the surface 



where L is the lateral size of the system, the scaling function f(x) ~ x@ for i < 1, and f(x) — > 
constant for x 3> 1. The scaling behaviors are characterized by the roughness exponent a, and 
the exponent (3. a features the long-time limit of the self-affine fractal structure as represented 
by the surface width w(L, t) ~ L a . At an early stage of growth, the surface width w (L, t) scales 
as t", where (3 = — and z is the dynamic exponent. 

The discreteness of a lattice is an indispensable ingredient in either experiments or computer 
simulations. As is known, due to this effect the crystalline 2D surface undergoes a roughening 
transition at a finite temperature. Such transitions of various systems have been observed 
experimentally in plastic and metal crystals. The theoretical understanding of this behavior 
is based on the Discrete Gaussian Solid-on-Solid model (DGSOS) |],|3||§] , in which the surface 
elastic energy is assumed to arise from the interactions between nearest neighbor atoms in a 
crystal. In addition, the periodic potential is introduced to account for the lattice discreteness. 

One-dimensional systems with only short-range interactions should not have a phase transi- 
tion. As a result, the DGSOS model in 1+1 dimension does not exhibit a roughening transition 
0. However, in nature there exist various kinds of long-range interactions between atoms, 
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such as the Van-der Waal force or the effective force due to bulk strain ||. The general form 
of interactions is found to scale as some power law with the distance between atoms. 

In view of this possibility, we consider a particular type of surface energy due to the long- 
rang interactions between atoms in a crystal @J^]. This model could be realized in vicinal 
surfaces whose steps are flat along one fixed direction. Here, we shall investigate the static and 
dynamic properties of a crystalline interface which undergoes the roughening transition in one 
dimension. Although qualitatively similar to the KT transition of the 2+1 dimensional surface, 
the 1+1 dimensional case exhibits different critical behaviors. 

The 2D DGSOS model with bulk or substrate disorder has been investigated recently |8|-[10|. 



The disorder was found to have non-trivial effects on both the dynamical and static proper- 
ties of the crystalline surface. The renormalization group (RG) analysis predicts a so-called 
super-roughening transition. The low temperature phase is dominated by disorder, and displays 
glassy behavior, identical to the vortex glass phase found in 2-dimensional dirty type-II super- 



conductors p"T|-|T4|j. Here we shall incorporate the same type of disorder into the ID interface 
with the elastic energy due to long-range interactions. 

The motivation of this study is twofold. First, we would like to investigate the behavior 
of an interface with the elastic energy arising from long-range interactions in two different 
environments- -pure and disordered. Secondly, current computer simulations of an 2D interface 
subject to quenched disorder are limited to modest spatial and temporal scales, and the present 
model could be used to probe more effectively the hydrodynamic limit. 

The elastic energy of the interface caused by the long-range interactions ||[7) considered 
here is: 
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h(x) — h(x')i2 K f°° dp 



H =^- dx dx ^ = - / ^\p\h(p)h(-p), (2) 
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where h(x) is the height of the interface at position x, and p is the strength characterizing 
the elastic energy due to long-range interactions between atoms. The temperature T has been 
absorbed into the Hamiltonian, and K = 2np/T. Powercounting shows that surface energy due 
to short-range interactions / dx(^) 2 is irrelevant. Therefore we will not take it into account 
here. 

The periodic potential inherited from the discrete nature of lattice is given by: 

a f°° 

H p — — cos [7/1(2;)], (3) 

CL J — 00 

where - is the coefficient of the leading term due to the discreteness, and 7 = x with a 
and b being lattice constants along and perpendicular to the base line, respectively. With this 
potential, the height of a surface should be regarded as a continuous variable within the interval 
(—00, 00). As in the DGSOS model, the higher harmonics are irrelevant near the critical point 
and can therefore be neglected. 

The action in Eq. (|2]) can also describe a variety of systems with quantum dissipations. A 
particular case is the electronic system involving a quantum wire with one impurity scattering 
center The Euclidean action of the electron at the impurity center is equivalent to the 



sum of Eq. (g) and Eq. (§). 

The dynamics of interface growth considered here is governed by [|3l,f4l,|l6f : 

dh(x,t) 5TL . . 

m =fJ.F- f j,— + r ) (x,t), (4) 

where t is time, F is the driving force, p, is the mobility, and 7i is the sum of 7io and 7i p . rj is the 
thermal noise which obeys (rj(xi, ti)r/(x 2 , £2)) = 2}J>8(xi— x 2 )S(t 1 — 1 2 ), where the coefficient 2p is 



chosen such that the system will reach the equilibrium states weighted by the Boltzmann factor 
e~( Ho+H p\ Consequently the fluctuation-dissipation theorem JlTj] will hold here in contrast to 



the case of far- from equilibrium surface growth [18], which cannot be formalized by Eq. 



To facilitate the analysis of Eq. (W) we employ the formalism of Martin, Siggia, and Rose 



(MSR) ]T9,20|. In addition to h(x,t) } the auxiliary field h(x,t) is introduced as well as their 
conjugate source J(x,t) and J(x,t). Upon averaging over the thermal noise the associated 
generating functional is obtained: 

Z[J,J] = J VhVhe Sca = J VhV~he So+Si »\ (5) 

where the effective MSR action S e s takes the form: 

/OO (^rp poo _ _ Q 

/ / dt{Jh + Jh + ^h{p,t)h{-p,t)-h{-p,t)[-h{p,t)-K^\p\h{p)}}, (6) 
-oo Z7T J-oo Ot 



oo roo 



S- mi = — / dx I dth(x,t) sm[yh(x,t)}. (7) 

d J— oo J —oo 

Here Eq. (§) will be treated by the RG technique. To perform the RG calculation, we follow 



the schemes of Amit et al. [21] and Neudecker [1G]. Regularization of the Feynman integrals is 



established by introducing the infrared cutoff M and ultraviolet cutoff a into the free correlation 

function C {x,t) = (h{x,t)h{0,0)) = JZo S^bITM e " /i( ^ H+M)|i|e " irawitn V 2 = x * + ° 2 - 
In the absence of the periodic potential, the 1+1 dimensional interface width scales as 



K l a v L 

where L is the linear size of the system, and the scaling function ln/(x) = — In 2 + Ei(47rx) 
with E\{-y) = -f™du^. 

I. Roughening transition in 1+1 dimensional systems: 
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To analyze Eq. @ we calculate the perturbative correction to several vertex functions up to 

2 

the lowest order in r = 1 — -^-^ and the second order in a, and use the minimal subtraction and 

Z7TA 7 



the standard prescription to obtain the flow equations of various parameters |i7|j23[| (details of 
the calculation will be published elsewhere). 

Under rescaling factors x — > e l x and t — > te' z (7 = lny where y is the rescaling factor), the 
RG recursion relations of various parameters are given by: 

^ = (* - 1 - (12) 

We first discuss static properties. Eq. @ holds true for all orders of a, since there are no 
conterterms proportional to \p\ in any order. Note that under the RG transformation only the 
terms proportional to powers of p, i.e. local in space, can be generated. Non-renormalization of 



K is also justified by using the real space RG calculation [T5|]. On the other hand, in the case of 
2D crystalline interface the surface tension does suffer renormalization. As shown in Eq. ( |TT]) 
the sign of r = 1 — = 1 — determines the fate of coupling constant a, where T c = b 2 p. 
When r < or T > T c , the flow of a under the RG transformation was driven to zero, and thus 
lattice effect can be neglected. In the long-wavelength limit, the surface width w(L,t = oo) 
scales as (lnL) 1//2 , as implied from Eq. (H). On the other side of the transition(T < T c ), as 
argued in Ref. 7, under the RG iteration coupling constant a flows to the infinity. On a 
sufficiently large length scale, a mass term generated by the cosine term will dominate over the 
surface energy strength K. The interface is pinned at the minima of the periodic potential, 



and is led to a flat phase as also occurs in 2+1 dimensional systems. Although there is no fixed 
point of a for r > 0, a(l) = a(0)e Tl should hold as long as I <C 1/r = Id- 

Next we turn to the dynamical properties of the interface in both phases. To facilitate the 
discussion, we give the expression of \x on the length scale / by integrating Eq. ([LI]) and Eq. (0): 

^(l) = ^(o) e 1} . (13) 

The linear response mobility (as F — > 0) can be obtained by letting I = oo in Eq. flUf). We 
find: 

7r 2 q,(0) 2 

Hm = (i(0)e * (14) 

for r < and hm = for r > 0. Here the parameters followed by (0) represent the correspond- 
ing bare values. In contrast to the roughening transition in the 2+1 dimension case, the linear 
response mobility decreases to zero as the temperature approaches the transition point from 
above. The higher order correction will not change this behavior qualitatively. In the case that 
the applied force F is finite, the movement of an interface is characterized by the non-linear 
response mobility. Since F is a relevant field, on some sufficiently large length scale I* ~ — In F, 
it will grow significantly. The pinning effect will be wiped out by the motion of an interface. 
Beyond this scale, the effect of the strength a, on the average, can be neglected. The mobility \i 
will scale normally, and its canonical dimension is zero. Consequently, the non-linear response 
mobility is yielded by stopping the RG iteration up to this scale: 

Mr~-lnF)~M0)exp[--^], r>0 (15) 

where C is a constant depending on bare values of other parameters. Meanwhile I* shoould be 
less than i.e. — InF < K Such a form of nonlinear response has been of interest to statics 
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of the vortex glass phase of type II superconductors. Here we see it can arise even in pure 
systems. 

In the region far below the transition point the results presented here in general will not 
hold. However, one can resort to the activated-dynamics mechanism. 
II. Depinning transition in 1+1 dimension: 

Crystalline surface growth with bulk disorder or two-dimensional substrate disorder has 
been addressed recently. Here we shall study the dynamics of a 1+1 dimensional interface with 
elastic energy in Eq. (0) subjected to substrate disorder. 

In the presence of substrate disorder, the periodic potential || changes to: 

a f 

Hdis = — J dx cos(j[h(x) + d(x)]), (16) 

where d(x) is the quenched fluctuation of the substrate height uniformly distributed within 
the lattice spacing b. The correlation of d(x) under question is presumed to be of short-range 
within lattice spacing a. In terms of the phase-like variable 6(x) = 7<i(x), we assume: 

(e i9{x) ) = 0, (17) 
(jO{x) e -ie(x')j =a s( x - x '). (18) 

As can be seen in Eq.(17), substrate disorder acts to destroy the structure of a regular lattice, 
thereby eliminating the roughening transition. However, as we show below, the correlation of 
disorder will lead to a new term which is relevant at lower temperatures. 

After averaging over disorder 9, the corresponding MSR action becomes: 

5i n t = f ^ x f ^ f dt'h(x,t)h(x,t)cos[yh(x,t)— r yh(x,t)], (19) 



oo J-oo J—oo 



where g(0) = a |^ , and again higher harmonics are discarded. 



Now we apply the RG analysis to the MSR action combining Eq. (^j) and Eq. ([19]) in the 
same spirit as we did the previous case. For the present theory the expansion parameters are 
t — 1 — -3^- and g. We derive the recursion relations for x — > xe l , and t — > te lz as following: 

dK , . 

IT = (20) 

^ = ( z - 1 - (23) 

As inferred from Eq. (p2|) above, the transition point is T pc = b 2 p/2, i.e. r = 1 — -f- < 0. 
The coupling constant g pertaining to disorder flows to zero under the RG transformation. 
Below T pc , g flows to a non-Gaussian fixed point g* = y. 

As with the previous case, we first discuss the static properties. In the high temperature 
phase, disorder is irrelevant and the surface width (or static height-height correlation function) 
is w(L,t = oo) ~ (InL) 1 / 2 . In contrast to the random DGSOS model, the marginal time- 
persistent parts p4j-[26[1 such as J J J dpdtdt' \p\h(p,t)h(—p,t') are not generated, since the term 



proportional to \p\ is spatially nonlocal. The form of the static correlation function will rely 
on the existence of such a term, as in the 2D case. Near the transition point the height-height 
correlation functions of two phases will not be of any difference in contrast to the 2D random 
DGSOS model. 

Next we turn to the dynamical properties. Above T pc , z = 1, and the linear response 
mobility remains finite and decreases to zero as the temperature approaches the transition 
point. The integration of the RG recursion relations gives the expression of the mobility at the 
scale /: 



M0=M0)[l + ^(e /i -l)]^ /2 . (24) 



and yielding the linear response mobility: 



k'r /2 

= /*o [2j7(Q)]7r/2 , as r ^ 0+. (25) 



Below T pc , due to the existence of an infrared stable fixed point g* , the dynamic exponent 
z assumes a non-trivial value 1 + |(1 — ^). The linear response mobility vanishes below T pc . 
Along the same line we derived Eq. ( |T5D the nonlinear response mobility below T pc is given 

by: 

H ~ MO)^ 1 "' 1 (26) 

to leading order in r . 

The physical implications of the RG analysis in the low-temperature phase are the following: 
Unlike the case of the 2D interface, the time persistent term in the ID case is not generated 
under the RG transformation. Physically the interface with long-range interactions cannot 
afford costly elastic energy to adjust itself to follow the short-range correlated substrate disorder. 
In the low temperature phase of the 2+1 random DGSOS model, the surface does follow 
disorder. In the present case, disorder is also relevant in the low temperature phase. In the 
presence of a force, the traveling interface on the average is attracive to the local minima of 
the random periodic potential. This implies a larger dynamic exponent and a smaller mobility 
in low-temperature phase than in the other phase. 

It is straightforward to extend the calculations for the long-range elastic energy with other 
powers in the denominator of Eq. @. The found behaviors are generic to systems with long- 
range interactions. The qualitative difference between systems with short- and long-range 
interactions is attributed to whether or not the interactions are renormalized. 
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In conclusion, we have investigated the dynamics of a 1 + 1 dimensional crystalline interface 



with long-range interactions. The governing equation is cast in terms of the MSR formalism 
and is analyzed by using the RG technique. In the presence of a regular lattice, above T c , 



the interface width is w(L,t) ~ ln(L/(-|)) ^ , and the linear-response mobility is /im = 



7j- a(u) 

fi(0)e 27 . Below T c the interface is flat. The nonlinear mobility is /i ~ /x(0) exp 



c 
7F^ 



In the presence of substrate disorder, above T pc the interface width is w(L,t) 



]n(Lf(j^j 1/2 . Below T pc , the surface width is w(L,t) ~ ln^L/(^ 



1/2 



with z = 1 + 



|(1 — y~) . The expression of the nonlinear mobility is also obtained for disordered case. 
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